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Abstract 

In this paper we discuss the effect of focal inhomogeneities on the 
globaf expansion of neariy FLRW universes, in a perturbative setting. We 
derive a generic iinearized averaging operation for metric perturbations 
from basic assumptions, and we explicify the issue of gauge invariance. 
We derive a gauge invariant expression for the back-reaction of density 
inhomogeneities on the global expansion of perturbed FLRW spacetimes, 
in terms of observable quantities, and we calculate the effect quantita- 
tively. Since we do not adopt a comoving gauge, our result incorporates 
the back-reaction on the metric due to scalar velocity and vorticity per- 
turbations. The results are compared with the results by other authors 
in this field. 

1 Introduction 

An essential difficulty which occurs when deahng with realistic cosmological 
models, is related to the fact that although the universe seems to be very close 
to FLRW at length scales of the order of the Hubble radius, the metric and 
matter content of the universe appears to be highly inhomogeneous at smaller 
scales. Since the realistic universe, with all its details at length scales small 
compared to the Hubble radius, is too complicated to handle in most calcu- 
lations, it seems desirable to extract those physical quantities which describe 
the large scale structure of the universe. However, when one tries to define an 
averaging operation for metrics, a number of difhculties occur. One of these 
difficulties is related to the fact that the Einstein equations are inherently non- 
linear, which makes it a nontrivial question to see how the Einstein equations 
constrain the dynamics of an averaged metric. Another fundamental problem 
which occurs when one tries to average metrics, is related to the fact that there 
is generally no direct physical significance in an averaged metric. Although 
this problem is usually ignored in the literature on averaging, it needs to be 
addressed before one can extend the discussion on averaging beyond an intu- 
itive level of understanding. The usual approach to averaging (see e.g., (l) - 
seems to be that one defines an averaging method, which is chosen on the basis 
of mathematical elegance or an intuitive notion of smoothness, and then one 



defines averaged physical quantities by means of the averaging operation which 
one has chosen. The objection against this approach is that if one calculates 
e.g. the averaged expansion of a perturbed FLRW universe, one can obtain vir- 
tually any result, by choosing an averaging operation which yields this specific 
result. In section ^ we explicify this problem, and we derive a generic linearized 
averaging operation for metrics which satisfies the condition that unperturbed 
FLRW is a stable fixed point of the averaging operation. It is shown that this 
generic linearized averaging operation for metrics can be expressed in terms 
of the spatial average of the perturbation of the spatial volume and goo in co- 
ordinates which are synchronous in the background. In section ^ we discuss 
the gauge problem and the choice of the background spacetime. The avera ged 
constraint equations are explicitly evaluated in section ^ and in subsection 4.1 



we derive an expression for the correction to averaged expansion due to density 
perturbations, in terms of the power spectrum of the matter. In subsections 
|4.2| and 4.3 we discuss the back-reaction on the metric due to matter velocity 
perturbations, and we show that vorticity perturbations may be important in 
the large wavelength limit. In section 4.4 we calculate the different corrections 
to the averaged expansion quantitatively by means of the observational data, 
and we compare our results with the results derived in previous works. In this 
paper, we adopt the convention that greek indices run from to 3, while latin 
indices run from 1 to 3. The metric signature is ( — h H — h), and the velocity of 
light, c, is set equal to one. 



2 The spatial average 

In this section we will consider the generic linearized averaging operation for 
metrics, for which unperturbed FLRW is a stable fixed point, and we show 
that this averaging operation has a universal limit when applied iteratively 
to perturbed FLRW spacetimes. We determine this limit explicitly, and by 
using the symmetry of the background, we show that the averaging of the 10 
components of the metric perturbation Sgpa- , can be expressed in terms of the 
spatial average of goo and in coordinates which are synchronous in the 

background. ^From now on the background FLRW spacetime will be called 
S, while the inhomogeneous spacetime is called S. Furthermore, we assume 
that S is coordinatized such that t represents the time coordinate which labels 
the hypersurfaces of homogeneity S in S, and E is coordinatized by x'' where 
i € {1,2,3}. We call a metric g^i/ or a metric perturbation Sg^i, spatially 
homogeneous and isotropic when there exists at least one coordinate system in 
which the components of g^j/ or Sg^i, are spatially homogeneous and invariant 
under spatial rotations. 

Let us consider the most general averaging operation A, which is a functional 
of metric perturbations about some background solution S, 

A6gf,^{x) = Tf,^{Sgpa{x)). (1) 

We will require the condition that unperturbed FLRW, in a gauge where the 
metric perturbation Sg^^, is spatially homogeneous and isotropic, is a stable 
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fixed point of the averaging operation A. This condition states that the av- 
eraging operation increases the spatial symmetry of the spacetime on which 
it works, assuming that this spacetime is sufficiently 'close' to FLRW, and it 
defines what we mean by averaging in this paper. 
It follows directly from this assumption that 

^^.(0) = 0, (2) 

for all X, since a nonzero value at the right-hand side of equation (j^) implies 
that unperturbed FLRW with the the same geometry as S, in a gauge where 
Sgfiv — for all is not a fixed point of A, which contradicts our assumption. 
The linear approximation to the averaging operation (|l|) , is given by 

A^^^^9^,,{x)^ fd^x' {x,x')6gp,{x'), (3) 
Js 

where the bi-tensor density fl^^(x,x') is defined as the functional derivative of 
J^fj^i, with respect to Sgpa-, evaluated at the point with coordinates x' in the 
background, i.e., 

(4) 

Sgp„=0,p=x' 

and we used condition (|^) which states that the zeroth order contribution in 
the expansion of A vanishes. 

The condition that unperturbed FLRW is a stable fixed point of the aver- 
aging operation A implies that the limit 

lim (5) 

n—^oc 

exists, and the quantity A°°Sgfj,i, must be spatially homogeneous and isotropic 
(we used the notation A^^^" to denote the n-times repeated operation of A'^^^). 

Note that the averaging operation A has two aspects; first it changes the 
geometry of the spacetime on which it works, and second it specifies a corre- 
spondence between points in the spacetime S, the averaged spacetime AS, and 
the background S. 

When one only requires that unperturbed FLRW is a stable fixed point of 
A, one constrains the way in which A changes the geometry of the spacetime 
on which it works, but one does not constrain the correspondence between 
points in the spacetimes S, AS and S. We constrain this freedom by impos- 
ing the stronger requirement that unperturbed FLRW, in a gauge where the 
metric perturbation Sgfj,i, is spatially homogeneous and isotropic, is a stable 
fixed point of A. This condition enforces that A does not generate 'pure gauge' 
perturbations when operating on unperturbed FLRW. 

Starting from equation (|^), and using the symmetries of the background 
spacetime S, it is shown in appendix B that the averaging operation A°° can 
be defined in terms of a spatial averaging operation which is universal, i.e., 

A^6g^,{t,x') ^{5gp,){t), (6) 



dgpaip) 
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where 

(%.>W= / rfVax/^ (7) 
Js(t) 



nP{x')n''{x') nf,{x)n,{x) + -hP" {x')h^,,{x) 



where nP denotes the future directed unit vector normal to S, and hP'^ :~ 
gBptr _j_ j^Pj^cr |g ^-|^g projection operator on S, and a denotes the distribution 
which is constant on E, and for which the integral over S equals one. Note that 
nPn"5gpa equals the perturbation of 500 in coordinates which are synchronous 
in the background (i.e., coordinates for which g^^ = —5°), while hP'^Sgpa equals 
the perturbation of the spatial volume element on E, to first order. It follows 
from this observation that the linearized averaging operation for metrics (0) , is 
effectively a spatial averaging operation for scalars, applied to Sgoo and 6g\ in 
coordinates which are synchronous in the background. 

An explicit realization of the spatial averaging operation for a scalar q{x), 
in the case where S is open, is given by 



{q{x)) = hm {q{x)){e) (8) 
lim N-\x,£) [ d^x' q{x')e{i ~ As{x,x')), 



where N{x,e) := x' yg^ 6(i — As(a;,a;')), and As(a;,x') is a distance 
measure between points x and x', ^ is a parameter with the dimension of length, 
and Q{x) = 1(0) for x > 0{x < 0). In the case where E is closed, (q) is defined 
analogously to expression (||), with N{x,i) =volume (E). 

It is shown in appendix A that the spatial average of a scalar function is in- 
variant under spatial gauge transformations, to arbitrary order in the expansion 
parameter of the gauge transformation. 

Notice that the spatial average (||) is only well defined when we make the 
assumption that perturbations q{x) are sufficiently small, such that the limit 
£ ^ (X) in equation (^) exists. It should be stressed that this assumption is 
nontrivial, and it is not automatically satisfied in general cosmological situa- 
tions, where perturbations are not necessarily bounded in amplitude and length 
scale. Indeed, since the observable part of our universe is restricted to our past 
light cone, there is no observational basis for the assumption that our universe 
is 'close' to FLRW at arbitrary large length scales. The usual way to deal 
with this situation is that one adopts a priori philosophical assumptions, such 
as the Copernican principle, to choose between different models which satisfy 
the observational data (see e.g., |^). Throughout this paper, we will adopt a 
version of the Copernican principle by assuming that perturbations are small 
enough such that the limit £ ^ 00 in equation M) exists. 



3 The gauge problem 

As is pointed out by Futamase in [Q, the observed matter density contrasts 
are of the order of unity at dimensionless length scales k of the order of 10~^, 
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where k denotes the fraction of typical size of the density fluctuation and the 
Hubble radius ru c/Hq. A rough estimation of the order of magnitude 
of the associated Newtonian gravitational potential, which we call e from now 
on, can be obtained by using the Poisson equation. For density contrasts of 
the order of unity we find e ~ k^, which implies a Newtonian potential 4> of 
the order of 10~^, suggesting that a perturbative approach might be adequate. 
At length scales of the order of the Hubble radius, the observable part of the 
universe appears to be highly homogeneous and isotropic, which motivates our 
choice for the FLRW metric as a background metric. 

Let us first briefly discuss some details concerning the spherical harmonic 
decomposition of perturbations about a background FLRW spacetime. 

The FLRW background metric can be written in the form, 

ds^ = dx^' dx" = (f) (- rfP + 77,^ dx' dx^ ) , (9) 

where rjij is the metric tensor for a homogeneous and isotropic three-space with 
curvature k, and f is a conformally scaled time parameter. We define the metric 
perturbation /i^^ by, 

fU..^fC + h^. , .g^"- - 3^^"^ - (10) 

and since g^^'g^^ = S^^, we have /i'"' = g^^Pg^^'-^hpa, and h^'^ = g^^'^'hp^, to 
first order. 

Copying Bardeen's notation in we define scalar, vector and tensor spher- 
ical harmonics Qn"* , Q^^] and Q|f ]j , respectively, which satisfy the Helmholtz 
equations (5,^''''''^^ + hnQn"^ — where p e {0, 1, 2} and | denotes the covariant 

derivative with respect to gfj. The vector harmonics Q^^^ are divergenceless, 

(2) 

while the tensor harmonics Q]j are divergenceless, symmetric, and traceless. 
We define traceless symmetric scalar harmonics Q^^jj by 

Q% - + and Q^:% - {kl - 3k)g(°) = 0, (11) 

and traceless symmetric vector harmonics Q^^lj ^-''^ defined by 

Q%--=-l^JQ% + Q^n]\d and gif -(fc2_2k)Q«=0. (12) 

The spherical harmonics are labeled by the parameter n € 0, Z+ (k G R^) 
in the case where S is closed (open). It is useful to define the hypersurface 
integration operation for scalars q(x) by 

m := {q (g(3)^/5(3¥), (13) 

which differs from the spatial average (^ by the volume element which is eval- 
uated in the background. As we show in appendix 2, the spatial average (|^) of 
a physical quantity is invariant under spatial gauge transformations, while the 
hypersurface integral ( p^ ) is generally gauge dependent at second and higher 
order in the expansion parameter of the gauge transformation. 
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The spherical harmonics Qn^ , (5^\^ and Q|f are orthogonal with respect to 
the hypersurface integration operation, i.e., 

{{Q^°^QfJ))= {{QMT))= {(0^,0^)'')) -5^n'. (14) 

and 

^((C^i'^^')) = 2((Q«.Q«^^)) =<5„„,. (15) 

Notice that the spherical harmonics are only to zeroth order orthogonal with 
respect to the spatial averaging operation (|^) due to a generally nonvanishing 
first order term which arises from the expansion of the volume element ^/g^ = 

V5^(l + /l + 0(/l2)). 

The most general representation of a symmetric 4x4 tensor /i^j, in terms 
of the complete basis of spherical harmonics is given by 



,o. = -a^E[B(o)Qi°i+i?WQW] (16) 



(2)n(2)i 

nijl 



where the coefficients An , Bn^ , Bn^ , H^p^ , H^p^ and H^^ are generally depen- 
dent on the conformal time parameter t. Let be the four-velocity associated 
with the frame in which the energy flux of the matter vanishes, then the three- 
velocity /vP associated with u'^, can be expanded as 

u^/u'^Y.'^vl?Q^:i+v^^^Q^::]], (17) 

n 

where qI^^ •= ~^n^Qi*ji' ^'^d u° = l/a(^ to first order, due to the normaliza- 
tion u^u^ = — 1. 

A gauge transformation is defined as a change in the correspondence be- 
tween points p in S*, and points p in S. The most general first order gauge 
transformation is the result of the coordinate transformation 

t-t + ^T„Ql")(x'^), (18) 



E 



(L(°)Q(")Xa:^)+iWgi^)X^'')), (19) 



in 5*, while the coordinates in S are fixed, and the correspondence between 
points with the same coordinates in S and in 5 is kept fixed. The coefficients T„ 
and Ln in expression ( |l8| ) and (^) are arbitrary functions of the conformal time 
coordinate i. Notice that T„ generates a change in the correspondence of the 
time coordinates in S and >§, while L„ generates a change in the correspondence 
between the spatial hypersurface coordinates on S and E. The changes in the 
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amplitudes of the metric tensor are calculated in the case of scalar perturbations 

= An ^ Tn T^, (20) 

a 

B^^B^+lJ^+k^^, (21) 

H^lI - H^l - (fcn/3)4°^ - -Tn. (22) 

<!=<i + fcn4°\ (23) 



where a dot denotes conformal time differentiation. For vector perturbations 
we find, 

r(i) - r(i) + r(i) h'^^^ - H^^^ + k L'-^^ (24:) 

while for tensor perturbations the trivial relation = H^n holds for all n. 
The matter velocity perturbation coefficients Vn'' and Vn^ , with respect to the 
coordinate frame, transform as, 

(25) 

where i € {1, 2}. Apart from the gauge freedom which is related to the mapping 
between points in 5 and 5, there is a gauge freedom related to the choice of 
the background scale factor a{t). A first order change in the choice of the 
background scale factor, 

h{t) = a{t) + D{t) (26) 

affects the spatially homogeneous mode of the trace part of the spatial metric 
by a change, 

rr(0) _ tt{0) , D 

^LO - ^LO + — ' (27) 

where is the coefficient which multiplies the spatially homogeneous trace 
mode in the expansion of the metric (p^), and 



(28) 



to first order. 

The approach in this paper will be based on a specification of the temporal 
part of the gauge (i.e., the correspondence between the time coordinates t 
in S and t in S), while maintaining covariance with respect to spatial gauge 
transformations (i.e., the correspondence between the spatial coordinates in 
S and in S). 

We stress that a fully gauge covariant approach is preferred to an approach 
which is based on a (partially) fixed gauge, since explicitly gauge dependent 
results generally point out nonphysical features of the calculation, including 
calculational mistakes. However, the intricateness of a fully gauge covariant 
calculation at second order makes such a calculation cumbersome (see e.g. ]lO|), 
and we will therefore follow an approach where the temporal part of the gauge 
is fixed, while maintaining spatial gauge covariance. 
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The temporal inhomogeneous part of the gauge is specified by imposing 
conditions on the coefficients in the expansion of the metric (|l6|) . The extrinsic 
curvature tensor of the constant-t hypersurfaces in S is given by, 



K' = 

3 



5] 



H, 



(0) 



(1) 



kB, 



(1) 



_L ff(2)n(2)« (29) 



By requiring that the coefficients An , Bn^ and in the expansion of the 

metric (^6|) satisfy the condition 



H-(O) 



-Ar. 



= 



(30) 



for all n, we specify a gauge in which the hypersurfaces of constant time t in S 
have spatially constant volume expansion K = ia/a? , as is clear by contracting 
expression (^). Condition (|30| ) specifies more or less uniquely a collection of 
spatial hypersurfaces in S (see ||ll[), but uniqueness is not required in the 
calculation which follows, since, as we will show in the following, our result for 
the average expansion of an inhomogeneous universe docs not in relevant order 
depend on the choice of the inhomogeneous temporal part of the gauge. 

Note that condition ( |30| ) does not constrain the choice of the time coordinate 
in S, and the correspondence between the time coordinates t in 5 and t in S. 
We specify the time parameter t in S*, up to the freedom of adding a constant, 
by imposing the requirement that the homogeneous component of A vanishes, 
i.e., 

^0 = 0, (31) 

for all times t. 

The choice of gauge ( |3T| ) implies that the background time interval coincides 
with the averaged proper time interval in 5, as measured by observers which 
are comoving with the spatial coordinates. 

The gauge condition ( |3l| ) can always be satisfied by performing a first or- 
der homogeneous gauge transformation. In order to clarify this statement, let 
us consider how equation (|3l| ) is affected by a homogeneous temporal gauge 
transformation. According to expression (po|), a homogeneous temporal gauge 
transformation with T — To, induces a first order change in the metric pertur- 
bation coefhcient Aq, 

io = - To - -To. (32) 
a 



The gauge condition (31) is satisfied by performing a gauge transformation of 
the form (pO), where 



dr a(T)Ao{T), 



(33) 



e 1 

and c is a constant of integration. The gauge condition (|l]) therefore deter- 
mines the homogeneous temporal part of the gauge, up to a constant of integra- 
tion c. According to the transformation law (22), the constant of integration 
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c in expression ( p3| ) affects tlie spatially homogeneous trace part of the spatial 
metric. This gauge freedom can be fixed by requiring that the homogeneous 
trace perturbation of the spatial metric vanishes, i.e., 

= 0, (34) 

for one time tc and for a fixed choice of the background scale factor a{t) at 
t = t^. 

Although we have completely specified the homogeneous temporal part of 
the gauge by imposing the gauge conditions and (jsj), the homogeneous 
trace perturbation of the spatial metric may still differ from zero for times t ^ 
ic- These perturbations are related to the freedom of choosing the background 
scale factor a{t) for times t 7^ ^c, as is clear from equation (p6|). When we 
require that condition (Q) holds at all times t, then it follows from equation 
( p6| ) and that the choice of the background scale factor a(t) is fixed for all 
times f € R. 

Recall that in section (j^) we derived the generic linearized averaging oper- 
ation for which unperturbed FLRW is a stable fixed point. It was shown that 
this linearized averaging operation which works on the ten components of the 
metric tensor, reduces to evaluating the spatial average of Sg\ and Sgoo- By im- 
posing the gauge conditions (^) and (|3^) , we specified a choice of background 
geometry by requiring that the spatial averages of Sg'^ and Sgoo both vanish. 
For this choice of gauge, the averaged spacetime equals the background space- 
time, and the averaging problem reduces to solving the averaged constraint 
equations for the background scale factor a{t). 

An explicit expression for the background scale factor a{t) in the gauge fixed 
by condition ( p^ ) is obtained by substituting the expression for the background 
metric and expression (jl^) for the perturbed metric, into expression (p8|). 
To first order we find, 

a'it) = ((^ )), (35) 

where g^'^-' = det{gij) and n — det{riij). 

Recall that condition (|30| ) fixes the inhomogeneous temporal part of the 
gauge, and the collection of spatial hypersurfaces on which the spatial average 
is evaluated. Since physical results must be gauge invariant to relevant order, 
one may question whether the freedom of choosing a family of hypersurfaces 
affects the result for the scale factor (|35|). It follows from the orthonormality 
relation (p^ and the transformation property (^2|) that the background scale 
factor ( p5[ ) is invariant to first order under inhomogeneous temporal gauge 
transformations. However, at order inhomogeneous metric perturbations do 
contribute to the background scale factor (|3^), and the gauge invariance of 
the scale factor a{t) therefore breaks down at order e^. Consistent with this 
limitation we will neglect terms of order in our calculation, while retaining 
terms of order e and j . 

Summarizing the content of this subsection, we completely specified the 
temporal and the spatially homogeneous part of the gauge, and the choice of 
the background, by imposing the gauge conditions ( |30| ) ( ^ and (|3^) on the 
metric coefficients A^^Bn^ and 
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4 Averaging the constraint equations 



The classical constraint equations on a hypersurface E are given by 

^(3) _^^2 _ K^^K'^ = 167rGp + 2A, (36) 

Kl.^ - K,, - 8^GJi, (37) 

where ; denotes the covariant derivative with respect to g^, and i?*-'^^ is the 
Ricci scalar associated with the induced metric gij , and 

p = T^'^n^n, , J, ^ -Tf^'^h.^n,, (38) 

where denotes the future directed unit vector normal to S, and :— 

In the constant-iiT gauge, defined by conditions (^, (|3l| ) and (|3^), the 
constraint equation ( |3^ ) takes the form 

where K'^^ := X*-' — hg^-' K is the traceless part of the extrinsic curvature tensor. 



In principle one could solve the constraint equation (39) for the time depen- 
dence of the scale factor a(t), while taking into account all linear and higher 
order contributions to the right-hand side of equation (|3^). However, this ap- 
proach is unnecessarily complicated, since all terms which do not have constant 
values on E must cancel on the right-hand side of equation (p9|) , since the left- 



hand side of equation (39) is constant on E. For the sake of calculational 
convenience, we will take the spatial average at the right-hand side of the con- 
straint equation (^9|) , without changing any physical aspects of the constraint 
equation; 

In order to solve equation (^ ) for the scale factor a(t), we need to evaluate 
the spatial average of the 3-curvature R'^'^\ and the energy density p, and the 
square of the traceless part of the extrinsic curvature tensor KijK^^ . We will 
calculate these quantities in the following three subsections. 

4.1 The averaged spatial curvature 

The spatial curvature perturbation SR^^^ can be expanded in terms of the 3- 
metric perturbation hij (see e.g. |9|), 

Ri3} = ^+SR'-'\ (41) 

where 

^6R+ 5^R + 0{h^), (42) 



and 



^^=<-^'l;> (43) 
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S^R = --h'^h^". +1/1^7;, l« - i/i/ij' +td., (44) 

where td. stands for terms which are total derivatives. Let us now evaluate the 
contributions to the averaged curvature perturbation (i?'^') for scalar, vector, 
and tensor modes in the expansion of hij . 

4.1.1 Scalar perturbations 

It follows from expression (^2|) that the lowest order contribution to the spa- 
tial curvature perturbation is given by g^^^SRij, which is order e/k^, but the 
spatial average of this contribution vanishes to order e/re^, due to the orthogo- 
nality relations (p^. The linear curvature perturbation g^^^SRij does however 
contribute to the averaged 3-curvature perturbation by a term of order e^/n^, 
i.e., 

{6R) ^ i| ^(fc^ - 3k)i/l°j(i/i':l + I H^l) + (45) 

n 

where we made use of the expansion of the volume element y/g^ = \/ 5^(1 -I- 
h + 0{h'^)), and the definition of the spatial average 

The quadratic term S^R in the expansion of the 3-curvature perturbation 
( ^2| ) contributes to the averaged 3-curvature perturbation by a term 

(^'^)^^ =-iE(fcn-3k) (10i/ir~^4r +lH'^yfl) + 0{e'/.% 

n 

(46) 

where we used the computer algebra package MAPLE to derive this expres- 
sion. Combining expression (|5|) and (^), we find an expression for the scalar 
contribution to the spatial curvature perturbation, 

(5i^(3))(o) ^ ^ ^(fc2 _ 3k)^2 ^ o{eyn% (47) 



where 



<^ + (48) 

is the gauge invariant amplitude which measures the distortion of the intrinsic 
geometry of the constant-X hypersurfaces. Using the expansion ( p2| ) for the 
spatial curvature perturbation, and the definition (^) of (j)hn, one finds that 
(j}hn is related to the first order spatial curvature perturbation by, 

5R^^Y.^kl-3k)cl,hnQi°^. (49) 



By substituting expression ([49| ) into the constraint equation (p6[), we obtain a 
simple expression for (f)hn in terms of the first order energy perturbation, 

(bhn = Jj^^^Gp e,„ + 0(eV«;'), (50) 

for all n, where e/i„ is defined as the density contrast in the constant- ii' gauge, 

e/m = Sp{kn)/p, (51) 
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and p denotes the background energy density. 

We would like to express the scalar contribution to the averaged curvature 
perturbation in terms of observable quantities. Since the averaged curvature 
perturbation ( p7[ ) is quadratic in (jj^n , we may use the constraint equation ( pO| ) 
to first order to determine (phn in terms of the fractional energy perturbation 
e/in- We obtain, 

(0) = 32{naGp)^ ^ + Oi^V^'), (52) 

where the sum (or integral when E is open) is taken over all possible values n. 
Expression (^2|) takes an especially simple form when expressed in terms of the 
power spectrum -P(fc), which allows the representation 



P,(fc) = ^^5(fc-|g„|), (53) 

where the subscript h refers to the constant- if gauge (see ag. or for 
more on power spectra). Combining expressions (52) and (|53| ) yields, 

(<5i?(3))(o) = 32(7r2Gp)2 + Oie^e^K^), (54) 

where 

J2 := ina^ / dkPh{k) (55) 



is an observable quantity known as the second moment of the power spectrum, 
and by absorbing a factor in the definition of J2 we restored physical units 
of length square. 



4.1.2 Vector perturbations 

Using the definition of the vector harmonics (p^, and the orthogonality rela- 
tions (p^, we find that vector perturbations do not contribute to the spatial 
curvature perturbation ( p2[ ) . This result may be expected, since it follows from 
expression ( p^ that one can always choose a gauge in which there are no vector 
perturbations of the spatial metric, and the vector contribution to the averaged 
spatial curvature perturbation (^) must therefore vanish in any gauge, due to 
gauge invariance of the averaged spatial curvature perturbation. 



4.1.3 Tensor perturbations 

Using equations ( [1^ ) , and expression ( ^2| ) for the second order expansion of the 
spatial curvature, it follows immediately that, 

n 

while the tensor contribution to the term (KijK^^) in the averaged constraint 
equation ( p3S|) fo llows immediately from the expression for the extrinsic curva- 
ture tensor (|29|). Although the tensor contribution to the averaged constraint 

— . . . (2) 

equations is easily calculated in terms of the coefficients the magnitude 

of this term has not yet been determined quantitatively by the observation of 
gravitational waves. 
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4.2 Averaged energy density 

In this subsection we will calculate the averaged energy density (p). In order 
to calculate the lowest order nontrivial contribution to the averaged energy 
perturbation, we will adopt the assumption in this subsection that the matter 
in the universe at late times after decoupling can be effectively described by 
the energy momentum tensor density for a pressureless perfect fluid, i.e.. 



(57) 



where is the four-velocity of the fluid, and pq is the energy density in the 
rest-frame of the fluid. The equations of motion for the fluid read. 



which implies 



where we used that V,, 



d, 



'9)~^dny/—9i and 



(58) 
(59) 

for a pressureless 



fluid. By using the spatial gauge freedom (|19|) we may set i?„ 



(0) 



^(1) 



0, 



such that \/—g — \/—goov g^^\ and the equation of motion (Mb takes the form 



(60) 



while in this gauge /vP equals the matter 3-velocity with respect to the 
normals to the constant-if hypersurfaces. The velocity four-vector u'^ can 
be written in the form 



1 



(61) 



where v\ := gijU^u^ equals to flrst order the square of the velocity three-vector 
/u^, and we used that u^u^ = —1. By substituting expression ( |6l| ) into the 
equation of motion (160), we find 



d_ 
di 



'g^po 



d 



= 0, 



(62) 



to flrst order. Using equation (62) and the definition of the spatial average (| 
we obtain 

hm ^ {{l + lyl)p,)(i) = 
£-^oo at I 



■£in(^lniV(x,£)) iS\^lvl)p,m 



lim iV~i(£) 



d 



dx' — y3^Vg(3)pou' e{e - Asix, x')), (63) 



where N{£) denotes the dimensionless quotient of N{£), and a constant with 
the dimension of a 3-volume. The second term on the right-hand side of equa- 
tion (^) vanishes due to Gaufi's theorem. Combining the remaining terms in 
equation (§3) yields. 



d 

lim —In ((1 

e^oo dt ^ 



A)po m 



d 

lim —In N{xA). 
i^oc dt 



(64) 
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By integrating equation (pih, it follows that 



where we used the gauge condition ( |34| ) . Formula (|65|) shows that the rest- frame 
energy density pQ, when integrated over a spatial volume element on T,{t) which 
is comoving with the matter flow, is not conserved for a pressureless fluid, while 
po{l + ^vf^) is conserved to first order. 

The spatial average of the energy perturbation Sp is obtained by expanding 
equation ( |57| ) for to first order, where we use (|6l|) and the gauge condition 
@. We find 

{p){t)^{il+vl)po){t), (66) 
which combines with equation ro5h. 



{P)it) ^ Pit) + ilpvDit), (67) 

to first order, where we used that (po)(^o) equals p{ta) when perturbations 
vanish at time to- Indeed, the lowest order contribution to the averaged energy 
density (|6|) , is given by the sum of the averaged rest-mass of the fluid, and the 
(nonrelativistic) kinetic energy of the fluid. Since vf^ is of the order of e, the 
lowest order correction to the averaged energy perturbation is typically small in 
the observed universe, but nevertheless signiflcant in the sense of the ambiguity 
which is related to the freedom of choosing a gauge and an averaging operation 
(see section ^) . 

It is interesting to note that there exists a simple relation due to Irvine 
and Layzer (see e.g., [Q) which relates W := 2TrGpJ2, where J2 is defined by 
equation (|5q), and the energy due to the peculiar velocity L :— {\pv\). For a 
pressureless fluid and nonrelativistic motions, it can be shown that -^{aW — 
aL) = La, which, assuming that the universe departs from small values of J2 
and L and relaxes to a nearly time independent bound state at late times, 
implies the Newtonian virial theorem L = W/2. 

Note that our result differs from a result derived by Futamase (see Q and 
[|j), where one finds a peculiar velocity contribution to the averaged energy 
density which is exactly twice as large as our result (|6^). This result seems to 
be based on the erroneous assumption that the integral of rest-frame energy 
density over a spacelike hypersurface is time independent (this is only true in 
a gauge where — vanishes) . In this case, equation (|5^) yields an averaged 
energy perturbation which is twice the result (|6^). However, this result violates 
continuity of the scale factor at the right-hand side of equation (|4^) when rest- 
mass is instantaneously and homogeneously converted into kinetic energy or 
vice versa. 

4.3 The squared shear contribution 

In this subsection we will evaluate the contribution of the term 

{K^,K^^), (68) 
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in the averaged constraint equations (40), for scalar and vector perturbations. 
The scalar and vector part of Kij are coupled to the matter current by the 
constraint equation (|37|), which takes the form, 

= SnGJ,, (69) 

when evaluated in the constant-ii' gauge. The matter current Ji is defined by 
expression (Bsh, and can be expanded as 



= iP + P)T. Q~ + ^i" ] , (70) 



to first order, where w^*^^ and denote the scalar and vector components of the 
velocity three- vector of the matter with respect to the normals to the constant- 
K hypersurfaces, Q^!^^ := ^k~^Q''^^^, and p and P denote the background energy 
and pressure density. 

By substituting the traceless part of the extrinsic curvature tensor (p9| ) , and 
the expansion ( |70| ) for Ji, into the constraint equation (|6^), we obtain 



- (fc„ - 3k/fc„) [4*^2 - A;„i3("'] = aG{p + P) v'^i^l (71) 



for scalar perturbations, and 



\ (fc„ - 2k/fc„) [ijW - A:„i3(i)] = aG{-p + P) v\2 (72) 



for vector perturbations. Expressions ([7l|) and (|7^) yield expressions for the 
scalar and vector traceless part of the extrinsic curvature tensor ( p9[ ) , in terms 
of the matter velocity, which can be used to evaluate the scalar and vector 
contribution to expression (p8[). For scalar perturbations we find. 



and for vector perturbations, 

ik^k^-) =2a-G\-p + Pf E (fc„_'£/fc„). - (74) 

The coupling between the matter current and the shear of the normals to the 
constant-X hypersurfaces, can be interpreted as the 'frame dragging' effect 
which occurs in the presence of moving matter (e.g., as in the region around a 
rotating black hole). It follows from expressions (|7^) and (^4|), taking into ac- 
count the normalizations of the scalar and vector modes (see expression (p^), 
that the matter current and Kij couple with different strength for scalar and 
vector perturbations. Furthermore, the strength of the coupling vanishes pro- 
portional to when fc„ — > oo. Since vj^^^ = 0(e), when velocity perturbations 
are generated by density perturbations at late times, it follows that expressions 
( [73| ) and ( |7^ ) contribute to the averaged constraint equations (^) by a term 
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of order en^ , which is neghgible compared to the leading order kinetic energy 
contribution discussed in section (4.2) for perturbations at length scales much 
smaller than the Hubble radius. 

However, for perturbations at arbitrary large length scales, the strength of 
the coupling grows proportional to S^-^ when S [0, where 6 := k^ — 3k for scalar 
perturbations and 5 := — 2k for vector perturbations. Note that since fc„ 
must be real for bounded solutions, the limit S I does not exist when k < 0, 
and the limit S I does not exist when k > since A;„ takes only discrete values 
in this case. 

Note that the divergent coupling between the metric and the matter velocity 
for S I and k = 0, is unrelated to the dynamics of the matter and metric at 
small scales and late times, since perturbations for which 5 <C 1 are typically 
larger than the Hubble radius, and must have a primordial origin. 

A natural question which arises is whether the divergence in equations ( [73| ) 
and @ for (5 i can be purely attributed to a large warping of the constant- 
K hypersurfaces, which can be removed by choosing another gauge. Indeed, 
it follows from expressions (29) and ( pl| ) that the scalar part of Kj can be 

set equal to zero, by a temporal gauge transformation with T ~ fc~^[ij^°^ — 

kBn''], but according to expressions (^), ( [2^ ) and (p3|), the intrinsic spatial 
curvature does diverge when 5 J, in this gauge. Furthermore, due to expression 
(24), the vector part of /Cj is gauge invariant, and the divergence in equation 
(74) is therefore independent of the choice of time-slicing. From the point 
of view of the matter, the most natural choice of gauge is a comoving time- 
orthogonal gauge, which is defined by the condition that the spatial coordinates 
are comoving with the normals to the constant-^ hypersurfaces (i.e., B^^^ = 
B^^^ = 0), and the scalar part of the matter velocity with respect to the normals 
to the constant-t hypersurfaces vanishes (i.e., -yC) — = — 0). According 
to expressions (^l|) and (|2^) , a gauge transformation from the constant- if gauge 
to a comoving time-orthogonal gauge is generated by T = In. this 

gauge, the scalar part of the shear of the matter coincides with the scalar part 
of ifj . By transforming equation (|7^) from the constant-X gauge to a comoving 

gauge, we find that the infra-red divergence of the scalar part of ifj has the same 
strength in both gauges, and its presence is therefore related to the presence 
of shearing matter. At first sight, a divergence of the shear of the matter for 
(5 1 seems to be inconsistent with the smallness of the velocity perturbations 
which are the source of the metric perturbations. There is no real inconsistency 
however, since the matter velocity perturbation is gauge dependent, and it 
might therefore be anomalously small in the constant-if gauge, without being 
in confiict with large matter shear perturbations. These observations show that 
the divergence in equations ( [73| ) and ( |7^ ) is of a physical nature. 

The absence of FLRW solutions of the constraint equations ( |69| ) when ho- 
mogeneous vector perturbations of the matter velocity are present, might seem 
peculiar, since solutions of the Einstein equations correspond to stable points 
of the action. At this point we should recall that we have limited our scope to 
FLRW background spacetimes, which are by definition spatially homogeneous 
and isotropic. In the presence of homogeneous matter velocity perturbations, 
our spacetime is no longer isotropic in the averaged sense, and there is no 
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FLRW background solution which is everywhere close to our perturbed space- 
time. A satisfactory description of homogeneous velocity perturbations about 
FLRW, requires the inclusion of background solutions which are homogeneous 
but not necessarily isotropic, and which include FLRW as a special case. These 
solutions are given by the Bianchi models of type V and Vllh, which include 
FLRW with k = — 1 as a special case, and type VIIq which includes FLRW 
with k = (see e.g., || - Q). 

4.4 The averaged expansion 

By substituting the expressions for the averaged curvature perturbation and 
the averaged energy density, which where derived in the previous subsections 



4.1 



4.2 and 4.3, into the averaged constraint equation (f40p, we obtain 



— = - — + -A+—{pvl) - {GpfJ2 75 

a* 3 3 6 3 

where J2 is defined by equation (pa), and the term g.w. denotes the contribution 



due to gravitational waves (see subsection [4.1.31) . We see that the averaged 
constraint equation ( fz^ ) takes the form of the standard Friedmann equation, 
plus a contribution due to the peculiar velocity of the matter, and a contribution 
due to the averaging of scalar and tensor metric perturbations. Let us now 
determine the magnitudes of the different contributions on the right-hand side 
of equation ( [zsj ) by means of the observational values for p and J2 . Estimates 
from the Lick and CfA catalogs ||l|| value J2 w 200/1"^ Mpc^, and p « 
1.88 X lQ~'^'^h?Q, g cm~'^, where ft, is a dimensionless factor which expresses the 
uncertainty in the value of the Hubble parameter Hq — 100ft km s~^ Mpc~^, 
and ft is believed to be between 0.5 and 0.85. Inserting these values in the 
different terms on the right-hand side of equation ([ZS]), one finds, 

1.14 X 10-35/^2^8-2, (76) 

3 

^(Gp)V2 ^ 1.0 X 10-39ft2r!2s-2, (77) 

and 

^(pi;2) ^1.3x 10-40/if]s-2, (78) 

where we used the relation v ~ (37rGpJ2)^^^ (see section ^^ ). According to 
equations (^- (|78|), and the constraint equation (|7|), the matter induced 
metric inhomogeneities act as a very small negative correction to the averaged 
energy density, equal to about 1.0 xVt part in 10"*, while the back-reaction 
due to the peculiar velocity of the matter acts as a positive correction to the 
averaged energy density, equal to about 1.2 parts in 10^. The small negative 
correction to the averaged energy density leads to a slight overestimation of the 
age of the universe io — f ^cT^ assuming that $7 = 1, equal to about 5 parts in 
10^. 
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5 Comparison with previous work 



The work on this paper started as a correction of the derivation by Futamase 
in iQ] on the points of the treatment of the gauge freedom (see section |[) 
and the choice of the averaging operation (see section ||). The paper was also 
inspired as an attempt to address the fundamental ambiguity which enters the 
calculation of any averaged metric through the freedom of choosing an averaging 
operation. 



In a recent independent paper by Russ et al \2l , the back- reaction due to 



density perturbations was calculated by using the relativistic Zel'dovich ap- 
proximation ||2^ in a comoving gauge. The expression derived by Russ et al 
for the back-reaction due to matter density perturbations agree in sign, but 
is roughly an order of magnitude larger than the result derived in this paper. 
Furthermore, a possible effect due to vorticity of the matter was ignored in 
that paper. It should be noted that direct comparison between the results by 
Russ et al and the results derived in this paper, is nontrivial due to the fact 
that the gauges used in the two papers are not related by a first order gauge 
transformation. Namely, a gauge transformation from the constant- if gauge 
to the comoving synchronous gauge requires i„ — ~Vhn due to equation (p5|), 
and Vhn — 0(e^/^) since t;^„ — 0(e) when velocity perturbations are generated 
by density perturbations at late times. By working in a constant-K gauge, 
we avoided the problem of a breakdown of the perturbative expansion which 
occurs in the comoving gauge (namely, since metric and matter density pertur- 
bations are of the same magnitude in a comoving gauge, metric perturbations 
get typically large at late times, even though the perturbations in the intrinsic 
geometry are generally small in the observed universe). 

Finally, we mention the paper by Buchert and Ehlers where one inte- 
grates the Raychaudhuri equation over a spatial hypersurface in a Newtonian 
background, and a globally vanishing correction to the averaged expansion was 
found. Although the Raychaudhuri equation is also valid in GR, the Newto- 
nian approximation enters the calculation where the correction to the averaged 
expansion is expressed in terms of a boundary term, which accounts for the 
difference between the Newtonian result and the nontrivial correction (fz^ ) to 
the averaged energy density derived in this paper. 



6 Conclusions 

We derived the generic linearized averaging operation for metrics starting from 
the requirement that unperturbed FLRW is a stable fixed point of the averaging 
operation. By a gauge invariant approach, we eliminated unphysical degrees of 
freedom in our problem, and we explicified the fundamental ambiguities which 
are related to the freedom of fitting the averaged spacetime to the inhomo- 
geneous spacetime. The leading order nontrivial corrections to the standard 
Friedmann equation are expressed in terms of the power spectrum of the mat- 
ter, and the effect is calculated quantitatively by means of the observational 
data. The dominant correction to the averaged expansion is caused by the back 
reaction of matter density perturbations, and leads to a slower expansion rate 
and an overestimation of the age of the universe by approximately 5 parts in 
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10^. The back-reaction of velocity perturbations, including vortical motion of 
the matter, appears to be negligible at small length scales. However, it was 
shown that the back-reaction of velocity perturbations can be significant in the 
large wavelength limit. 
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8 Appendix A 

In this appendix, we will discuss the relation between the volume element in 
the hypersurface integral, and gauge invariance at second and higher order in 
the expansion parameter of the gauge transformation. 

Let (f) he a, one parameter group of diffeomorphisms : R x E — > E, which 
is defined by the condition that (/)a=o is the identity, and the curves 4>x{p) are 
integral curves of a vector field ^ in S (see e.g. Q and for the mathematical 
details which are involved) . A gauge is specified by choosing a mapping between 
points p in S, and points p in E. Assuming that a choice of gauge has been 
made, then a one parameter group of gauge choices is obtained by mapping 
the points 4>\{p) in E to points p in E, for all A G R (the more generic case 
of a one parameter family of mappings of points in the background and the 
perturbed spacetime, is discussed in pO[ |, but there is no need to introduce this 
complication in the derivation which follows). 

Let us consider a scalar function q{x), which lives in E (such that its value 
in a point p in E is fixed, while its value in a point p in E depends on the 
choice of gauge). The spatial average and the hypersurface integral of q{x), are 
related by 

(9) = ((?(5^'V5^)^)), (79) 

where we used the definition (|l^). The integrand at the right-hand side of 
equation ( |79| ) is gauge dependent, and may be expanded in powers of A about 
A = 0, i.e., 

k — QO , ^ 

g(ff(^)/ff^)^(A,p-) = J2 If4 9(^(^75^)^ (80) 

fe=0 

where £| denotes the fc-th order Lie derivative with respect to ^, evaluated in 
p. By substituting the expansion ( pO[ ) in the integrand at the right-hand side 
of equation (^ , we obtain 

k—oo , 

(g)(A) - (g)(0) = E TT {{Clq{g^'^/g^)h}. (81) 

k=l 
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For k — 1, the contribution to the right-hand side of equation (gl]) is evaluated 
using 

c^q = Cq■,^ (82) 

and 

/:?(5('V5^)*=CM5^'V5^)^ (83) 

where ; denotes covariant differentiation with respect to gij . Combining equa- 
tions (||) and (11) yields, 

Ciq{g^'^/g^)^ = {qeU9^'^/9'')K (84) 

and 

((('Zf);.(5('V.9'')^»-0, (85) 

due to Gauss's theorem. The fc = 2 contribution to the right-hand side of 
equation ( |8l] ) is obtained by making the substitution q — s- {qC);i in expression 
(^4[), and for arbitrary k G Z+ the same result follows by induction. Since the 
terms at the right-hand side of equation ( pl| ) vanish for all fc, we established 
that the spatial average of a scalar function q is gauge invariant to arbitrary 
order in the expansion parameter A. Applying the same analysis as above to 
the hypersurface integral of a scalar field q{x), we find, 

((9»(A) -(('Z))(0)= fc! (86) 

which depends on A, due to equation (^2|), unless g is a constant on E. A similar 
derivation, where we reverse the roles of the spacetimes E and E, shows that the 
hypersurface integral of a scalar field q{x) which lives in E is gauge invariant, 
while the spatial average of q{x) is gauge invariant iff q{x) is constant in S. 
It follows from these observations that the spatial average of a perturbation 
6q := q{x) — q{x) is gauge invariant iff q(x) is constant on E. 

Note that Futamase in 0] uses the hypersurface integral as a spatial av- 
eraging operation in the calculation of second order effects, while he does not 
consistently fix a gauge in these papers (namely, he assumes a comoving syn- 
chronous gauge, and constant expansion on the hypersurfaces of constant time 
coordinate). 



9 Appendix B 

In this appendix we derive the decomposition of the generic linearized averaging 
operation 

:= lim (87) 

n — >oo 

in terms of the spatial average {5g^i,){t)^ which is uniquely defined. Note that 
the existence of the limit ( |87| ) implies that 

A'^'^hl. = (88) 

for arbitrary spatially homogeneous and isotropic perturbations (up to the 
freedom of diffeomorphisms acting at either side of equation (pq)). 



20 



Without loss of generahty, a spatiaUy homogeneous and isotropic perturba- 
tion Sg*^^ about S can be written in the form 

^9*^^ = (t>i{t)nf,nt, + (f>2{t)hf,i,, (89) 

where hf^i, := gj^^ + n^fi^, and denotes the timelike future directed vector 
in S which is orthogonal to E, and which is normalized with respect to the 
background metric g^^^, and 0i(i) and 02(f) are arbitrary functions of t. 

When we substitute expression ( |89| ) for (5(7*^ and expression (i) for 
into condition (^8|), we obtain, 

= (j)i{t)nf_,{x)ni,{x) + <i)2{t)hi_ii,{x), (90) 

for arbitrary functions 0i {t) and 02 (i) . Equation ( pO| ) holds for arbitrary 0i [t) 
and 02 (i) iff 

/ d^x np{x')na{x') f^l{x,x') = 5{t' - t)np{x)n^{x), (91) 

and 

/ rfV V(x') /;^(x,x') = 5{t' - t)hp,{x). (92) 

Expression (^) shows that fp'^{x,x') is proportional to a delta distribution 
5{t — t'). It follows from this observation that A^^^ can be naturally defined in 
terms of a linearized spatial averaging operation A^s'\ i.e., 

=i(i)(t)%„ (93) 

where ^i^-* is defined by, 

ii'^%.= / d'^x' fPZ{t,x\x^') 5gp.{x% (94) 
■/set) 

and 

/;,-(^,x^x^'):- / dt' fPl{t,t',x\x^'), (95) 

and At' is chosen such that t E At' . At first sight, the decomposition of the 
linear averaging operation A^^^^ in terms of a spatial averaging operation which 
is defined on a collection of spatial hypersurfaces might be surprising, since the 
choice of a collection of spatial hypersurfaces Y,{t) in S is gauge dependent. 
It will be shown in section (||) that although the choice of E(t) in S is gauge 
dependent, the linearized spatial averaging operation (|9^) is to first order gauge 
independent. 

Assuming that the limit (H) exists, then by substituting expression ( |9^ ) into 
expression (g7|) one finds that the limit 

{6gp,) := hm (96) 
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exists. We will show that the limiting spatial averaging operation which is 
defined by equation (|9^) is universal. 

Expression (M) and the definition (Bq) imply that, 



(<55^,):=lim / d'x' f;rit,x\x'') dgp^ix'^), (97) 



where fj^^"^ is defined in terms of f^"^ by induction over n, 

/;r(^\ ^'^) = / d^qCpit, x'\ q^) x\ g'), (98) 

and fjj^'' := fP^. Let us now try to determine the limit 

■■= lim f^r- (99) 

An explicit calculation of /^'"^, using the definition definition ( |98| ) for fJ^P"' 
and starting with arbitrary realizations for /^J^, would be quite cumbersome, 
but fortunately it appears that the symmetries of the background spacetime S, 
and the stability condition (^ determine f^'"^ completely. 

Recall that we required that the limit ^ converges to a spatially homoge- 
neous and isotropic metric perturbation for arbitrary perturbations which 
implies that 

(%.>(a:0= / d^x' f^/-it,x\x'^)SgpAx'^)^Sg;,, (100) 



for all X, where we used expression (|9^) and Sg*^ has the form (|89|). If expression 
(100) holds for arbitrary perturbations 5gpa{x''^), it also holds for arbitrary 
perturbations Sgp„{x'^ + c'), where G R. By abs orbing the constants c* into 
the coordinates a:*, one finds that expression ( |lOQ| ) remains unchanged under 
the substitution 

f^P-{x\x'^)^ f^/-ix\x'^-d). (101) 

Furthermore, since the right-hand side of equation ( |100| ) is spatially homoge- 
neous by requirement, we find that the left-hand side of equation (100) must 
be also invariant under the substitution 

f^rix\x'^) -> f^rix' + d\ x"), (102) 



wher e G R is arbitrary. Since equation (100) is invariant under ( 101 ) and 
(102) for arbitrary perturbations Sg^i,, we conclude that ff^P'^ is (up to the 
freedom of performing diffeomorphisms) constant on S when regarded as a 
distributi on (i .e., neglecting sets of Lebesque measure zero). Furthermore, since 
equation ( |100|) holds for arbitrary Sg^i^, the distribution f^^'^ix,x') must be 
proportional to a tensor of the form ( p^ in the point x, thereby fixing the fiv 
dependent part of fj^P'^- We may therefore write 

f^nx,x') =g{''{x')n^{x)n,{x) + gf^'' {x')h^,{x), (103) 

where g^" {x') and gf^'^ {x') arc spatially homogeneous tensor densities in x' , and 
we used expression ( |89| ) for g*^. 
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A similar argument, using the invariance of equation ( |100| ) under the group 
of spatial rotations (using that T does not explicitly depend on x) , shows that 
the bi-tensor density f'^f"^{x,x') is isotropic with respect to the indices a and 
p, which implies that the tensor densities g'l'^ [x') and 52*^(2;') in expression 



^ nP{x')n''{x') +a2\/g^hP''{x') 



(104) 



(103) are of the form, 

gl (a; ) = ar 
and 

(a;') = a3\/^ nP{x')n''{x') + ai^/^hP" {x') (105) 

where 5^^^ denotes the real space volume element, which follows from requiring 
spatial gauge invariance at higher orders (see appendix A), and the factors a„ 
(n e {1,2,3,4}) are constant on E. Substituting expressions (104) and (105) 
in expression (|103|) yields 



f^n^.x) ^^){ainP{x')n''{x')n^{x)n,{x) + aJiP%x')h^,{x)), (106) 

where we used expressions ( |9l| ) and ( |9^ to show that the terms proportional 
to a2 and vanish. 



By substituting expression (106) for /^'"^ into condition (IOC), where we 
set 6gpfy equal to 5g*^^ d efine d by expression (|8^), wc find that the constants 
ai and 04 in expression (|106| ) must satisfy the condition 



d^x'^ 



E(t) 




1. 



(107) 



Expression (107) shows that the constants ai and 3a4 are equal to (volume(I]))~^ 
when S is closed, while in the case when E is open, ai and are defined in 
a distributional sense by condition (107), and by the condition that ai and 
are constant on E(t). 



By substituting expression (|10^) into expression ( |100| ) we obtain the explicit 
expression for the spatial average, 




l- 



[nP{x')n''{x') n^,{x)n,{x) + -hP%x')h^,,{x)] Sg, 



(108) 



(109) 



Note that nPn'^ 5gp„ equals the perturbation of goo in coordinates which are 
synchronous in the background (i.e., coordinates for which g^p = — 15"), while 
hP'^Sgpa equals the perturbation of the spatial volume element on E, to first 
order. 

Summarizing the derivation in this appendix, we showed that the general 
linearized averaging operation which is a functional of metric perturbations 
about FLRW, and for which unperturbed FLRW is a stable fixed point, has 
a unique limit when applied iteratively to perturbed FLRW. Furthermore, we 
showed that this linearized averaged operation is naturally defined in terms 
of a spatial averaging operation which works on goo and the spatial volume 
perturbation in coordinates which arc synchronous in the background. 
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